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Abstract 



m 

' We consider the spin calculations in arbitrary and diagonal spin bases 

(N ' 
O 

^ ■ 1 Introduction 

Qh! At the present time, the conventional method of direct calculation of the squared modulus 
r-| ! of amplitudes is practically not used in calculating of reactions. This is due to the fact 
that even an insignificant increase in the number of particles participating in a reaction 
as well as the account of their polarizations sharply complicate calculations. The idea of 
direct calculation of amplitudes followed by their quadrating is fairly obvious. However, 
such an approach entails a number of specific problems. Therefore, in the last few decades 
publications devoted to the development of particular aspect of this method have been 
appearing in the press (see references in reviews [0). 

These works can conventionally be classified according to the following directions: a 
covariant or an incovariant approach is used; calculations are carried out in an arbitrary 
or a particular spin basis. Among the first works devoted to this topic are the works of 
F.I. Fedorov and his scholars. For details see They are devoted primarily to covariant 
calculation of amplitudes in an arbitrary spin basis. For instance, in 0] a method that 
permits using in calculating amplitudes the main "trick" of the quadrating procedure was 
proposed. Namely, reduction of calculations to the computation of traces from 7-matrix 
operators. 



*Presented at the Vlth International Schoohseminar "Actual Problems of High Energy Physics". 
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Let us illustrate this approach by an example of a fermion "sandwich" which can be 
given in the following form [Q - arbitrary operator): 

u''\p\s')Qu''{p,s) = Ti Qu''ip,s)u'''{p',s') . (1) 

The main difficulties of this approach are the obtaining of an explicit form of the 
transition operator 

u''ip,s)u'''ip',s') , (2) 

as well as the solution of such problems as the calculation of the phase factors, the removal 
of singularities, the calculation of the exchange diagrams, and the account of the specificity 
of the spin configuration a' = —a, and so on. Some of these difficulties can be overcome 
by fixing a certain manner the spin basis. Transition to another basis is realized with the 
use of the Wigner D-functions. However, this apparatus is exceedingly cumbersome and 
incovariant. 

The present paper shows how, using amplitudes (|l|) in the basis in which they have 
the simplest form, to obtain expressions for the squared modulus of an amplitude in 
an arbitrary spin basis. The paper considers the following questions. In Section 2, one- 
particle states are considered. It is shown that consistent use of tetrads formalism permits 
operating with a complete set of bispinors as well as reducing the operations of matrix 
operators to tensor ones. 

In Section 3, matrices that make it possible, using the amplitudes calculated in one 
spin basis, to calculate the squared modulus of amplitudes in another spin basis. 

Section 4 considers the introduced into the diagonal spin basis (DSB). In the DSB, 
the bispinors of the initial and final states have a common set of spin operators and 
interrelated in the simplest way. As a result, the central problem of this approach - the 
obtaining of a simple and compact solution. The diagonal amplitudes have a clear physical 
interpretation and are the best "building bricks" for calculation of squared modulus of 
amplitudes in arbitrary spin basis. 

In Section 5, some examples of calculation in DSB are considered. 



2 One-particle states. Tetrads 

The one-particle state is defined by 4- vectors, momentum p and the axis of spin projections 
s. And ps = 0. This condition is fulfilled if s is determined in terms of the arbitrary vector 

q 

{pq)p — m^q 



{pqY — m?q^ 

Define 



(3) 



Sq = v = —,s^ = s, (4) 
m 
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then Sq = —si = 1, SqS^ = 0. Any vector having the form aso + (3s3 belongs to the 2-plane 
(50,53) in Minkowski space. The orthogonal 2-plane can be determined with the aid of 
two tensors 

gl'' = v^^v" - sf's" , = e'"'f"'vps^ , (5) 

where £0123 = 1- If the arbitrary vector r lies incompletelly in the 2-plane (s^Ss), then 
the vectors {g is Minkowski tensor) 



st={9f-9n^, s'^ = -efy: = ^r,{gf - ^M.)r. (6) 



satisfy the conditions = si = —1, S1S2 = and both vectors are orthogonal to sq and 
S3. Thus, the tetrad sq, S3, si, S2 form ortho normalised basis in the Minkowski space. In 
all spin reactions, vectors Si and S2 only occur in combination Si ± is2- From follows 

r 

si + ias2 = T„ — , T„= g\\- g- iae\\ . (7) 

It is clear that the tensor plays the special role in spin calculations. Besides the 
apparent properties T* = T^a, = = ~'2'T„, the relation 

rj-\fiurj-iaf3 rpflf^rpau f Q\ 

holds for it. From this relation follows 

Ti^^r^r^Tf = rlTf . (9) 
If in d^) we replace the vector r by r', then from (|^ and (|^) we obtain 

si' + 2as2' = e*"'^(si + 2as2) , e*-'^ = . (10) 

Thus, variations of r lead to rotation by angle (p in the 2-plane (si,S2), and in the spin 
relations an additional phase factor appears. Having fixed r, we thus fix the tetrad in the 
Minkowski space and thus can write in explicit form the whole set of spin operators^ : 

S3 = 7553, = ^(si +icrs2) , (11) 
which act on the bispinorQ u'^{p., s,r) in following way: 

T.su^'ip.s.r) = au^ip^s^r); ,^2) 
E_^u''{p,s,r) = u'''{p,s,r) , E^'u'^(p, s, r) = . 



-'^ 75 = i7°7^7^7^ can be given in form 75 — isiS2S3So. 

^ We will keep symbols s and especially r only where necessary. 
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The second relation represents the known phase agreement: the action of the spin flip 
operator does not lead to the appearance of an additional phase factor. Note the following 
important fact. Having fixed the tetrad, we obtain the possibility of writing for the 4- 
component bispinor four equations, namely, two equations for the spin flip operator are 
added to the Dirac equation and the equation for the spin projection. 

Let us consider of the consequences of this situation. We give the tensor g in the form 

gt^i^ = yM^^ _ s^'s'' - s^^s"; - . (13) 
Then from ([TT]), (|TB|) follows 



(14) 



= vi^v - 75 {s^^s + {s^i - ^crsi^)J:^ + (s^ + «as^)S_ J . 
Using the Dirac equation and ([12|) , we have 

^^u'^ip, s) = {v^ - as'^^^Kip, s) - (s^, + iasi^)^5U-^ip, s) . (15) 
Repeated action of the 7-matrix on ([T5|) gives 



[v ■ s]^'^75 + ■ s] j u'^i^p, s) 
- [(si + ias2) ■ {as + 75^)]'"'" u'^i^p, s) , 



(16) 



where the notation for tensor is introduced 

[a ■ b]"" = a^'b" - bf'a" . (17) 
Any operator Q can be expanded by complete set of Dirac matrices 

r = {l,75,7^757^c^^n • (18) 

Thus, formulas (p!5|), (|T6|) provide the possibility of reducing the action of an arbitrary 
matrix operator on the bispinor: Qu'^{p,s) to expansion by a complete set of bispinors 
ii^'^(p, s), 75ii^'^(p, s). The expansion coefficients will be tensors whose rank is equal to 
the number of free Lorentz-indices in the operator Q. 

3 The transit matrices between the different spin 
bases 

Let us assume that we know the set of amplitudes in a certain spin basis 

Ms'5 = u''{p',s')Qu\p,s), (19) 
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but the initial state is represented so that the axis of the spin projections is the vector /. 
Thus, we are interested in the amphtudes 

M's,^ = u''{p',s')Qu''{p,l), (20) 

or to be more exact, in the final analysis, the squares of their modules. 
Using the equation 

J2\{l + v){l + 5^,s)u"{p, I) = u^ip, I) (21) 

we obtain the expression for the squared module of the amplitude ( pO]) (arguments in 
bispinors are omitted): 

= ^ u^'Q^il + v)il + 6^^,s)u'^u'^^il + ^)(1 + ^275^)01^'' 

S1S2 \ ' 

5i(52 (5i(52 

where 

JCs.s, = u^'u^'ifu^^ = Ttu''u''u^'u^' . (23) 

Thus, if we know the explicit form of matrix (p3|), then using (^9|), we calculate the 
squared modules of amplitudes (|20|) . 
If Si = S2 = S, then 



ICss = 4 Tr(l + v){l + a75/)(l + v){l + dj,s) = ^ (1 - a6{ls)) . (24) 

Using the operators of spin fiipping 2'^5{si — iSs2), we obtain the expression for 
the configuration 61 = —62 = S: 

^s-s = 7^Tr(l + v){l + cr75/)75(si - ids2){l + v){l + S'y^s) 

1 ^2 . . a (25) 

= — Tr75(si - i6s2){a-^^l + absvl) = --(si - ^^52)/ • 
io / 

In deriving (^), we used the representation (53 = s, Sq = v) 

75 = isihhso . (26) 

/^From the explicit form (p4|), (|25| ) of the elements of the matrix /C^^^j it follows that it 
can be written with the aid of Pauli matrices af. 

^5rS2 = (1 - (^(^iilSi))s^s2 ■ (27) 
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(P7[) can be given in a covariant 4-dimensional form. To do this, we introduce the isotropic 



4- vector 

L = {l-a{ls{)) (28) 
and take into account that ctq = 1- Than 

K. = a^L^ . (29) 

It should be noted that 1 is nothing but the square of the 4-velocity of the particle, and 
/ and S3 are the axes of spin projections in different spin bases. 

Note that into (pT]), (^) enter only the vectors of the tetrad in which we calculated 
amplitudes (|l^). It can be shown that the finite expression ( p2|) is independent of not 
only the phase vector r, but also the vectors si and S2, since they enter into (^) in 
combinations SiS^ + S2S2 = S3S3 — gf^'^ and £^'^'"^SipS2o- = [vo ■ 83]^'^. 

Using formulas (|2^, (p5|), (|27D, ( |29| ) we can construct a matrix /C for any fermion 
participating in the reaction. And the representation of (|27|) , (^9]) provides the possibility 
of using, if necessary, the Fierz transformation. 

Note the following important fact. Formula (|2^) can also be used in the case where 
the beam is partially polarized and the states are not pure. To do this, it is enough to 
carry out substitution a/ — > a, |a^| < 1 in (|2^), (p^), ([281). a is the vector of partial 
polarization (see p). 

To investigate the polarization phenomena, besides the squared modules of ampli- 
tudes (l2^) , it is necessary to know the values of A^crA^l^J^ since various polarization 
characteristics are expressed in their real and imaginary parts. 

Thus, in the general case, it is necessary to know the matrix 

^^J^^ = Trw'^M'^M^^M''^ , (30) 



and at a2 = cri, we obtain (22) 



It can be shown that the matrix JC is of the form {5ai3 are Kroneker symbols) 

^sll' = \ {5.,..h,8, - ar^af'^ks,)] . (31) 

It should be recalled that with respect to 5i and 62 summation with the amplitudes 
calculated in the spin basis with the tetrade {v, S3, Si, S2} is made. 
In conclusion, note that the formula 

M'.^MZ, = E Ms,!^sll'Ml (32) 

<5i S2 

represent a new way of describing the transition from one spin basis to another without 
using the extremely cumbersome and noncovariant apparatus of D-Wigner functions. 

■^The projection of only that particle whose polarization we are investigating is given in the amplitude 
designations. 



6 



4 Diagonal spin basis (DSB) 



In any process, only an even number of fermions can participate. If the number of free 
fermions is 2n, each diagram contains n open fermion hues and is described by a structure 
consisting of contractings between n fermion "sandwiches" of the type 

^'''{p\s')Q^-{p,s). (33) 

For the fermion line $ are bispinors u (see (|1|)) and for the antifermion line - bispinors v\ 
for annigilation and for creation of an pair $ are bispinors u and v. 
For certainty, we shall consider the structures (|l]), (H). 

We shall not consider the group aspects of introduction of the DSB. We only note 
that in DSB the vectors s and s' are chosen such that they belong to the 2-plane {p,p') 
p p' 

or (v,v'); V = — , v' = — -. Satisfying this requirement, we obtain 
m m' 

{vv')v — v' , {vv')v' — V 
s = \ , s' = - \ . (34) 

^Jivv'y - 1 v(w)2 - 1 

Thus, the support vectors (see (^) for the on- and off- states of the fermion line are 
vectors q = v', q' = —v. With such a choice of signs in special reference systems vectors 
s and s' coincide with the direction of the 3-momentum of the initial particle and are 
opposite to the 3-momentum of the final particle. By special reference systems is meant 
the Breit system for fermion or antifermion lines (t-lines) and s.c.m. for the pair being 
annihilated or created (s-lines). 

/^From ( |5^ ) it can easily be seen that in the derived reference systems diagonality gets 
the meaning of helicity, with S = X, 6' = —A'. Thus, the DSB is covariant description of 
the helicity in the special reference systems. Exactly in these systems helicity the pairs 
of particles as well as such notions as non-flip and flip amplitudes have a clear physical 
meaning. Indeed, in helicity basis in arbitrary reference system the spin of the initial 
particle is projected on the 3-momentum p and that of the final particle on p', than what 
can be said about the non-flip or flip process? These notions in helicity basis have a only 
marking meaning. 

In our opinion, neglect of this fact is the main reason why the process of constructing 
operators (0) convenient for calculation has been extended to decades. Attempts to con- 
struct covariant operators (Q) in helicity basis look unreasonable. The helicity of massive 
particle is a "bad" quantum number, since it is not invariant at a Lorentz transformation. 
Any declaration of covariancy of operators (|^) in the helicity basis usually is a disguised 
transition to the special reference system. It should be noted that other fermion pair 
"droop", since each pair has its own special reference system. The only exception is 
e~^e~ — > type reactions. No wonder that just these reactions are chosen, as a rule, 

for examples of calculations of processes. 
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Before to go to the construction of tetrads for the initial and final states, we introduce 
into the 2-plane {v,v') two v- and v'- symmetrized, orthonormalized vectors 

v + v' v-v' 
^0 = ^ , = ^ ; (35) 



W ±1 

V± = \ . (36) 



From (PD, (^ it follows that 



(7y = v^v — s^s = V V — s s = UqUq — njjrig , 

/ Vj^ V— 

iFiom. (^, ( pTf ) it follows that the choice of common phase vector r' = r leads to the 
coincidence for both vectors of the tetrads lying in the orthogonal 2-plane, i.e. 

= s'i = s'^, = (gf - , = 5^ = ^'^ = -ef^ . (38) 

II II j.^ 

Next, let us coincidence the plane Lorentz transformation that transforms v to v'. In 
the representation of the group SL{2, C), it is of the form 

A{v - ^0 = , (39) 

and A^A+ = v'. From (03) it also follows that AsA+ = s'. Transformation (|39|) does not 



change vectors lying in the orthogonal 2-plane. Thus, the Lorentz transformation (p9|) 
converts the tetrad of the initial particle into the tetrad of final particle: 

A{v v') {v, s, n, 712} A^(f — > v') = jv', s', n', , (40) 

and this in turn means that the relation between the bispinors of the initial and final 
states in DSB is of the form 

u\p', s') = A{v v')u\p, s) . (41) 
In the DSB we choose the normalization 

u^{p)u^{p) = u^\p')u^' {p') = 1 . (42) 



Then the relation (^) describes the cases m! ^ m too. 
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To restore the generally accepted normalization, it is necessary to multiply the ampli- 
tudes calculated in the DSB by factor 



n J2mam[ , (43) 



where n is the number of open fermion lines. 

Using the Dirac equation and (p!4D, we can write relation (^) in different representa- 
tions 

u\p\ s') = "^^u^p, s) = hou\p, s) = {V+ - 6^,V^)u\p, s) . (44) 



Formula (|4^) makes it possible to express in the DSB the explicit form of the transition 
operator (^ in terms of the projective operators of the initial (or final) state 

u\p, s)u\p, s) = ^{v+l){l + d-f,s) . (45) 

As relation (@), the transition operators (H) can be given in different form. Some of 
them are shown below. 

Au'ip, s)u\p\ s') = + 1) - + 1) = 

= {V+ + VI) (1 - <575non3) + no + 575^^3 = (46) 

= + - {V+ + 5'^^V-){hQ + 575n3)^ (no + 575n3) , 

4.u\p, s)u-\p, s) = 

= — ^ + 1 -— r - / - ^ ^ - > / K + 1 = 

r_L \2V- vv' + 1 2V+ vv' — l J (-47^ 

= 75 (VV + 575 - no)(ni + i^ng) = 
= 75 (v+ + 575 V^_ - 2 ("-0 + '^75"-3)) (^1 + iSn2) . 

In deriving formula (|47|) from (^61) , the spin flip operator ( pT]) is used and the relation 
valid for an arbitrary orthonormalized tetrad also hold: 

no(ni + i6n2) = S'j^n^ini + i^ng) , (no + S'jr^n^Y = 2(1 - S-f^noU^) . (48) 

Since any interaction operator can be expanded by the complete set of Dirac matrices 
(p!8D, it is tempting to calculate the matrix elements of this set in DSB. From (|1]), (|i6|), 
it follows that 

u\p',s'){i-n,;r-n5r;<^n^^\p.s) =_ 
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u-\p', s') {1; 75; 757^^; s) = 

= {0 ; ; 5V4ni + 18712^ ] -V+in, + z^na)'^ ; 5 [^3 ■ {fii + t5n2)r} . ^ ^ 

The matrix elements (^), (0) can be interpreted as spin characteristics of exchange 
particles under the scalar, pseudoscalar, axial and tensor interaction, respectively. 

/^From ( ^OD it follows that the exchange (emission, absorption) of a scalar or pseu- 
doscalar particle does not change the fermion spin projection in DSB. Thus, the choice 
of DSB not only leads to simple calculation formulas, but also adequately reflects the 
physical meaning of the processes being described. To support this thesis, we consider in 
DSB the matrix elements of nucleon current 

■J'/ = u''ip', s') (^Fi(g2)7, - ^^^,uQ'^ ^\p^ s)-q=p'-p. (51) 

lYiom (H), (|0|) we have 

J'f = GeUo^ , Jj-' = SV-GMini + iSn^)^ . (52) 

It is noteworthy that in DSB formfactors Fi and F2 independently go in Saks form- 
factors having a clear physical meaning: 

Ge = Fi + -^F2 , Gm = Fi + F2, (53) 

Ge and Gm describe the distribution of the electric charge and the magnetic moment 
of nucleon, respectively. Thus, from ( p2D it follows that the situation where the non- 
flip processes are responsible for the electric interaction and the flip spin processes - for 
the magnetic interaction is covariantly described in DSB. It should be noted that the 
space-like {q"^ < 0) virtual photon has a scalar or a circular polarization, respectively. 

In helicity basis, grouping of (|53|) occurs only in Breit's system in which the Fourier 
transformations becomes three-dimensional since 

exp{—iqr) = expi{qf— qot) = exp{iqr) . 



Formulas (|^), ( ^Tj) described the fermion t-line. To describe the antifermion line as 
well as the s-line for the annihilating and creating pairs, one should make use by the 
relation 

v\p, s) = -Sj5U-\p, s) , (54) 

which relates the particle and antiparticle bispinors. 

If in (^), ( ^7D we restore by the recipe ( ^5D the normalization u^{p)u^{p) = 2m and 
perform the limiting transition m ^ or/and m' 0, we will obtain the transition 
operators for processes involving massless fermions found in 0. 
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For calculating concrete processes, it may be convenient to utilize a formalism, in which 
the basis spinor u^{nQ, n^; rii, 712) common for initial and final bispinors and satisfying the 
conditions 

hou\no, ris) = u\no, 71,3) , -f^hsu^no, 71,3) = 6u\no, 113) , (55) 

is used. ^ 

It is easy to seeF] that plane Lorentz transformations A(?7,o ^ v) = - — and 

V/2(VV + 1) 

A(?T,o v') = , ^ — — change the basis spinor tetrad to tetrads of the initial and 
V/2(VV + 1) 

final states respectively. Therefore, the transition operators (|46|), (^Tf) can be given in the 
form 

Au\p, s)u''ip',s') = ,^^(v+l)u\no,n,)u''in,,ns)iv' + l) = 

1 ^+ + ^ (56) 

In this equality the Bouchiat and Michel relation is used. 

In ([56|), vectors v and v' can be expanded in vectors no and 77,3 with the aid of relations 



In conclusion of this section, we give the recipe for calculating in DSB exchange di- 
agrams. As an example, we consider the electron-electron scattering. For certainty, let 
particles 1, 3 and 2, 4 are paired. Then the exchange diagram has the structure 

u''iP4Hu'^iPl)u''iP3h''u''iP2) = TTjy'ipi)u'^'ip3)ru''iP2)u''ip4) , (57) 

i.e. it is expressed in terms of the transition operators entering into direct diagram. 



5 Examples of calculations in DSB 

In Section 4, it is shown how the main task of the approach based on the use of formula ([l|) 
is solved in DSB in compact, covariant form. Namely, the problem of finding in explicit 
form the transition operators (Q). In each particular reaction, a further simplification of 
calculations can be achieved by simplifying the form of the interaction operators Q. Let us 
demonstrate this by an example of calculating the processes of radiation of bremsstrahlung 
photons. Let the initial fermion in the state u^{p, s, r) emits a photon with a momentum 
k and helicity A. We do not concentrize the spin projection vector of the fermion so far. 
Choose the photon polarization vector in the form P] 



1 



e^{k) = ^^{ei\+iXe\\- Xa)k, (58) 



"'For this the relation ~ — is used. 

V+ + 1 V- 
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eil f^'^ = vf's" - s^^v" , e\\ = ^e'"'f"'e\\^^ . (59) 

It is easy to see that vector (^) meets the necessary requirement for the polarization 
vector: ke'^ = (e^)^ = 0, e^e~^ = —1. /^From the last equation in (|T2]) and condition (p!0|) 
it follows 

i-yeiiku^ip, s) = -cr-yigw - g)ku''{p, s) , (60) 

therefore, takmer mto account that e v = 1= , we obtam 

y^k^ 

e^^ip, s) = e^v{l + \j5Kip, s) . (61) 

In amplitude, the operator is followed by operator p — k + m. Using ( pT| ) and the 
commutation relation pe = 2pe — ep, it can easily be seen that 

(p-k + m)e\''(p, s) = 2e^p (l - —(1 - A75)A;) u'^ip, s) . (62) 

V 2m ) 

The fact is rather surprising that if n photons with equal helicity A are emitted from 
the fermion line, then the structure avalanche- like increasing in complexity acquires in 
this approach a very trivial form 

(/n + m)en ■ ■ • (/i + riiYx u^^p, s) = 
= 2"/ne„ ■ ■ ■ /lei (1 - ^(1 - A75)A:) u'^ip, s) , ^^^^ 

here f i = p — ki — ... — ki is the momentum of fermion between the photons ki and /sj+i, 
k = ki -|- . . . -|- kfi. 

We now consider the following problem. In those diagrams where a photon is emitted 
from the final fermion line, we shall still use vector (p5|). However, we shall not obtain 



a structure similar to (|62D because tensors (|59|) were constructed for the initial fermion 
line. This problem disappears in DSB where these tensors coincide (P7|). 

If a photon is emitted from another fermion line, then by virtue of gauge invariance 
we can make the following substitution in these diagrams: 

(see 0) where the vector e' constructed by the recipe of ( ^8)) on basis vectors of another 
fermion line. 

If we pass to limit of massless fermions, we will obtain the known results CALCUL- 
group [0 , including the factorization of the sum of the contributions of all the diagrams 
obtained on permutation of photons in (|63|). 

If considering the radiation of gluons, it is necessary to take into account details 
associated with the gauge invariance in non-Abel theory [|1^ . 
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In the processes involving two bremsstrahlung photons, the case Ai — A2 — A is 
described by vector (|58D . For the polarization Ai = — A2 = A, it is possible to construct 
a polarization vector common for two photons by formulas analogous to (pSD- In this 
case, the amplitude structure takes on the form of a one-photon process and fact of the 



participation of two photons is only reflected by the kinematics factor |T3 



In [0 and review a number of processes going to field of an intense laser wave 
are considered. One of these processes, + 717* = + 7 (n is number of coherently 
absorbing laser photons) is basic in the operation of 67 and 77 colliders [0], the others 
are calibration and test ones. Calculation at level of amplitudes makes it possible to 
investigate the nonlinear and spin effects. And the spin effects are especially significant 
exactly in the nonlinear processes. 

If reaction involves W^, bosons, the structure of the transition operators (^6|) , (W^) 
is considerably simplified. This is due to the relations of the type 

(1 ± 75) (no + ^75^3) = (1 ± 75) (no + Shs) , 

i.e. these operators are expressed in terms isotropic tetrads, and the relation (1 ± 75) (1 =F 
75) = reduces to zero some of the terms. 

If two bosons participate in the process, their tetrads are plotted in the same manner 
as for a pair of fermions. In so doing, the circular polarization vectors of bosons coincide. 
With the participation of tree bosons, as, for example, in the reaction e~^e~ — ^ W~^W~ , 
it is convenient to choose for them a common support vector q = ki + k2 + = Pi + P2 
and a common vector 712- And 77-2'^ ~ s^^''"kiyk2pk^(j. This will considerably simplify the 
interaction operator Q whose structure comprises up to 5 7-matrices. 
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